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Abstract 

When  connected  in  small  amounts,  the  impact  of  distributed  generation  (DG)  on  distribution  system  stability  will  be  negligible.  However, 
if  its  penetration  level  becomes  higher,  distributed  generation  may  start  to  influence  the  dynamic  behavior  of  the  system  as  a  whole.  This 
paper  presents  a  mathematical  representation  of  a  solid  oxide  fuel  cell  plant  that  is  suitable  for  use  in  distribution  system  stability  studies. 
The  model  is  applied  to  a  distributed  utility  grid  that  uses  a  solid  oxide  fuel  cell  plant  as  distributed  resource.  Examinations  include  transient 
stability  and  voltage  stability  of  the  system. 
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1.  Introduction 

Distributed  generation  (DG)  is  electricity  generation  sited 
close  to  the  load  it  serves,  typically  in  the  same  building 
or  complex.  The  DG  embraces  a  palette  of  technologies  in 
varying  stages  of  availability,  from  entrenched  to  pilot.  It  is 
sometimes  called  a  “disruptive”  technology  because  of  its 
potential  to  upset  the  utility  industry’s  apple  cart. 

Most  likely,  fuel  cells  will  be  a  dominant  DGs  [1,2].  These 
DGs  are  dynamic  devices  and  when  connected  to  the  distri¬ 
bution  system  they  will  affect  its  dynamic  behavior.  Hence, 
several  researchers  are  working  to  develop  dynamic  models 
for  these  components  [3-10]. 

This  paper  develops  a  generic  dynamic  model  for  a 
grid-connected  SOFC  plant.  The  model  is  defined  by  a  small 
number  of  parameters  and  is  suitable  for  planning  studies. 

The  steady-state  power  generation  characteristics  of  the 
plant  are  derived  and  analyzed.  Understanding  the  transient 
behavior  of  SOFC  is  important  for  control  of  stationary  util¬ 
ity  generators  during  power  system  faults,  surges  and  switch¬ 
ings. 

Voltage  regulation  is  one  of  the  main  problems  in  the  dis¬ 
tribution  systems,  especially  at  the  much  far-end  load  and 
in  the  rural  areas.  Voltage  regulation  and  maintaining  the 
voltage  level  are  well-known  problems  in  the  radial  distri¬ 
bution  network.  Several  techniques  have  been  applied  by 
implementing  many  devices  in  the  distribution  network  to 


*  Corresponding  author.  Tel.:  +34-953-026518;  fax:  +34-953-026508. 
E-mail  address:  fjurado@ujaen.es  (F.  Jurado). 

0378-7753/$  -  see  front  matter  ©  2003  Elsevier  B.V.  All  rights  reserved, 
doi:  10. 1016/j.jpowsour.2003. 11.001 


solve  these  problems.  The  most  common  devices  and  tech¬ 
niques  used  are  transformer  equipped  by  load  tap  changer, 
supplementary  line  regulators  installed  on  distribution  feed¬ 
ers,  shunt  capacitor  switched  on  distribution  feeders  [11]  and 
shifting  transformers  towards  the  load  center  [12].  A  multi¬ 
ple  line  drop  compensation  voltage  regulation  method  that 
determines  tap  positions  of  under-load  tap  changer  trans¬ 
formers  is  proposed  in  [13]  to  maintain  the  customers’  volt¬ 
ages  within  the  permissible  limits. 

The  model  derived  is  based  on  the  main  equations.  It  is 
developed  in  the  Laplace  domain  and  transient  simulation 
is  done  using  a  software  developed  based  on  the  MATLAB 
package. 

The  paper  is  structured  as  follows.  Section  2  presents 
a  review  of  transient  stability.  Some  basic  concepts  of 
voltage  stability  are  introduced  in  Section  3.  Section  4  de¬ 
scribes  the  SOFC  model.  Section  5  briefly  discusses  the 
utility-connected  inverter  control.  Section  6  depicts  some 
simulation  results.  Finally,  conclusions  are  presented  in 
Section  7. 


2.  Transient  stability 

Transient  stability  is  a  term  applied  to  alternating  current 
(ac)  electric  power  systems,  denoting  a  condition  in  which 
the  various  synchronous  machines  of  the  system  remain  in 
synchronism,  or  in  step  each  other.  Conversely,  instability 
denotes  a  condition  involving  loss  of  synchronism,  or  falling 
out  of  step  [14]. 
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Nomenclature 

Fuel  cell 

Eo 

ideal  standard  potential 

F 

Faraday’s  constant 

Ifc 

fuel  cell  current 

Kh2 

valve  molar  constant  for  hydrogen 

Kr2o 

valve  molar  constant  for  water 

Ko2 

valve  molar  constant  for  oxygen 

Kr 

constant,  Kv  =  No/AF 

No 

number  of  cells  in  series  in  the  stack 

P 

real  power 

p* 

set  point  for  the  real  power 

Pi 

partial  pressure 

€ 

input  fuel  flow 

4 

fuel  flow  that  reacts 

r 

ohmic  loss 

m-o 

ratio  of  hydrogen  to  oxygen 

R 

universal  gas  constant 

T 

absolute  temperature 

Te 

electrical  response  time 

Tf 

fuel  processor  response  time 

u 

fuel  utilization 

Vic 

fuel  cell  voltage 

m2 

response  time  for  hydrogen  flow 

JH20 

response  time  for  water  flow 

JO2 

response  time  for  oxygen  flow 

Inverter 

E 

load  bus  voltage 

E* 

set  point  for  the  load  bus  voltage 

Lt 

inductance 

Q 

reactive  power 

Q* 

set  point  for  the  reactive  power 

V 

inverter  output  voltage  space  vector 

XT 

reactance  =  Ljco 

<$p 

angle  between  and  V'e 

s; 

angle  reference 

s; 

flux  vector  associated  with  V 

flux  vector  associated  with  E 

f* 

flux- vector  reference 

Distribution  network  stability 

Dxf 

system  Jacobian 

H 

inertia  constant 

Pc 

electrical  power  output,  in  p.u. 

Pi 

initial  load 

Pl 

active  load  at  a  bus 

Pm 

mechanical  power,  in  p.u. 

P  n  L 

load  distribution  factor 

Po 

active  power  base  load 

X 

state  variables 

APl,  A0l 

changes  in  load  demand 

5 

rotor  angle,  in  electrical  radian 

A 

changes  in  the  system  loading 

OOQ 

nominal  speed,  in  electrical  rad/s 

Hence,  for  a  simplified  intuitive  description  of  transient 
stability,  a  power  system  may  be  regarded  as  a  set  of  syn¬ 
chronous  machines  and  of  loads  interconnected  through  the 
transmission  network.  Under  normal  operating  condition,  all 
the  system  machines  run  at  the  synchronous  speed.  If  a  large 
disturbance  occurs  the  machines  start  swinging  with  respect 
to  each  other,  their  motion  being  governed  by  differential 
equations.  Depending  upon  the  power  system  modeling,  the 
number  of  such  first-order  differential  equations  is  lower 
bounded  by  twice  the  number  of  system  machines,  but  may 
be  orders  of  magnitude  larger. 

DG  units  normally  supply  power  to  the  local  load  centers 
but  the  excess  power  could  also  be  exported  to  the  regional 
power  grid,  adding  to  the  capacity  and  stability  of  the  grid 
system. 

The  key  to  interconnection  is  the  safety  of  the  people 
who  have  to  clear  faults  on  the  line,  and  protecting  the  DG 
generator  from  feeding  into  a  low-impedance  fault.  A  fault 
will  knock  DG  off  the  system,  requiring  it  to  be  resynchro¬ 
nized  with  the  grid.  There  are  various  means  to  enhance  the 
transient  stability  performance  of  the  distribution  systems 
[15-17]. 

Fast  valving  is  one  of  the  most  effective  and  economic 
means  of  improving  the  stability  of  a  power  system  under 
large  and  sudden  disturbances.  Fast  valving  schemes  involve 
rapid  closing  and  opening  of  thermal  turbine  valves  in  a  pre¬ 
scribed  manner  to  reduce  the  generator  acceleration  follow¬ 
ing  a  severe  fault  [18]. 

During  steady-state  operation  of  a  power  system,  there 
is  equilibrium  between  the  mechanical  input  power  of  each 
unit  and  the  sum  of  losses  and  electrical  power  output  of 
that  unit.  The  problem  arises  when  there  is  a  sudden  change 
in  the  electrical  power  output  due  to  a  severe  and  sudden 
disturbance.  The  severity  is  measured  by  drop  of  this  power 
to  a  very  low  or  to  zero  value  and  a  consequential  sud¬ 
den  acceleration  of  the  machines  governed  by  the  swing 
equation: 


2 H  d2(5 

A  A  —  ^>m  I  ’' 
co o  d tz 


where  8  is  the  rotor  angle,  in  electrical  radian,  Pm  the  me¬ 
chanical  power,  in  p.u.,  PQ  the  electrical  power  output,  in 
p.u.,  H  the  inertia  constant,  in  MW  s/MVA,  coq  the  nominal 
speed,  in  electrical  rad/s. 

From  (1),  it  is  apparent  that  the  decrease  in  the  mechanical 
power  has  the  same  impact  on  the  rotor  angle  swings  as  that 
of  increase  in  the  electrical  power  output.  Fast  valving  has 
a  function  of  reducing  the  mechanical  power  input  to  the 
turbine  and  so  the  generated  power. 

A  change  in  load  demand  (AFT,  AQl)  causes  corre¬ 
sponding  change  in  both  bus  voltage  magnitude  and  the 
power  output  of  the  generator.  The  voltage  magnitude  value 
is  calculated  in  the  power  flow  analysis.  The  change  in  the 
generator’s  active  power  output  results  in  system  frequency 
variations  through  the  swing  equations. 
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3.  Voltage  stability 

Voltage  instability  spans  a  range  in  time  from  a  fraction  of 
second  to  tens  of  minutes.  Time  response  charts  have  been 
used  to  describe  dynamic  phenomena  [19,20].  Many  elec¬ 
tric  system  components  and  controls  play  a  role  in  voltage 
stability.  Only  some,  however,  will  significantly  participate 
in  a  particular  incident  or  scenario.  The  system  character¬ 
istics  and  the  disturbance  will  determine  which  phenomena 
are  important. 

3.1.  Mid-term  voltage  stability 


The  time  frame  is  several  minutes,  typically  2-3  min. 
Operator  intervention  is  often  not  possible.  This  scenario 
involves  high  loads,  and  a  sudden  large  disturbance.  The 
system  is  transiently  stable  because  of  the  voltage  sensitiv¬ 
ity  of  loads.  The  disturbance  (loss  of  generators  in  a  load 
area  or  loss  of  major  distribution  lines)  causes  high  reactive 
power  losses  and  voltage  sags  in  load  areas. 

3.2.  Longer-term  voltage  stability 


The  instability  evolves  over  a  still  longer  time  period  and 
is  driven  by  a  very  large  load  buildup,  or  a  large  rapid  power 
transfer  increase.  The  load  buildup,  measured  in  MW/min, 
may  be  quite  rapid. 

Several  papers  have  shown  the  direct  relation  between 
saddle-node  bifurcations  and  voltage  collapse  problems,  e.g. 
[21,22].  Saddle-node  bifurcations,  also  known  as  turning 
points,  are  generic  codimension  one  local  bifurcations  of 
nonlinear  dynamical  systems  of  the  form 


dx 
d  t 


fix, 


where  x  e  Rn  are  the  state  variables,  X  e  R  is  a  partic¬ 
ular  scalar  parameter  that  drives  the  system  to  bifurcation 
in  a  quasi-static  manner,  and  /:  Rn  x  R  — >  Rn  is  a  non¬ 
linear  function.  System  exhibits  a  saddle-node  bifurcation 
at  the  equilibrium  point  (xo,  Ao),  he.,  fix o,  Ao)  =  0,  if  the 
corresponding  system  Jacobian  Dxf  |o  =  Dxf(x o,  Ao)  has 
a  unique  zero  eigenvalue,  and  some  particular  transversality 
conditions  hold  at  that  equilibrium  point. 

The  state  variables  x  and  nonlinear  function  /(.)  are 
typically  defined  in  terms  of  the  quasi- steady  state  phasor 
models  used  in  transient  stability  studies.  Thus,  generator 
angular  speed  deviations  and  phasor  bus  voltages,  magni¬ 
tudes  and  angles,  are  usually  an  important  part  of  vector  x. 
The  parameter  A  is  typically  used  to  represent  changes  in 
the  system  loading,  regardless  of  the  load  model  used. 

Therefore,  typically  one  starts  from  an  initial  stable  op¬ 
erating  point  and  increases  the  constant  power  loads  by  a 
factor  A,  until  the  singular  point  of  the  linearization  of  the 
power  flow  equations  is  reached.  The  loads  can  be  defined 
as: 


PL  =  P0il+XPnL) 


where  Po  is  the  active  power  base  load,  Pn  l  represents  the 
load  distribution  factor  and  Pl  the  active  load  at  a  bus  L  for 
the  current  operating  point. 


4.  Solid  oxide  fuel  cell 

There  are  several  types  of  fuel  cells  being  developed  for 
a  variety  of  applications  and  these  have  been  extensively 
discussed  in  the  open  literature.  Unlike  other  variants,  the 
SOFC  is  entirely  solid  state  with  no  liquid  components.  Op¬ 
eration  at  elevated  temperature  is  needed  to  achieve  the  nec¬ 
essary  level  of  conductivity  in  the  cell’s  solid  electrolyte 
for  it  to  operate  efficiently.  With  an  outlet  temperature  in 
the  range  of  900-1000  °C,  the  efficiency  of  the  cell  alone  is 
about  50%. 

Typically  the  fuel  cell  system  consists  of  SOFC  genera¬ 
tor  modules  in  a  parallel  flow  arrangement,  with  the  number 
of  standard  modules  being  determined  by  the  plant  power 
requirement.  The  SOFC  generator  module  embodies  a  num¬ 
ber  of  tubular  cells,  which  are  combined  to  form  cell  bundle 
rows,  several  of  which  are  arranged  side  by  side  to  make  up 
the  complete  assembly. 

The  high  efficiency  of  the  system  means  that  less  carbon 
dioxide  is  generated  than  in  contemporary  power  plants.  In 
addition,  the  fuel  is  oxidized  electrochemically  without  any 
interaction  with  atmospheric  nitrogen  so  negligible  amounts 
of  nitrogen  oxide  are  discharged  to  the  environment. 

This  paper  provides  a  basic  SOFC  power  section  dynamic 
model  used  for  performance  analysis  during  normal  opera¬ 
tion.  Some  control  strategies  of  the  fuel  cell  system,  response 
functions  of  fuel  processor  and  power  section  are  combined 
to  model  the  SOFC  power  generation  system. 

The  chemical  response  in  the  fuel  processor  is  usually 
slow.  It  is  associated  with  the  time  to  change  the  chemical 
reaction  parameters  after  a  change  in  the  flow  of  reactants. 
This  dynamic  response  function  is  modeled  as  a  first-order 
transfer  function  with  a  5  s  time  constant. 

The  electrical  response  time  in  the  fuel  cells  is  gener¬ 
ally  fast  and  mainly  associated  with  the  speed  at  which  the 
chemical  reaction  is  capable  of  restoring  the  charge  that  has 
been  drained  by  the  load.  This  dynamic  response  function 
is  also  modeled  as  a  first-order  transfer  function  but  with  a 
0.8  s  time  constant. 

With  aid  of  the  inverter,  the  fuel  cell  system  can  supply 
not  only  real  power  but  also  reactive  power.  Usually,  power 
factor  can  be  in  the  range  of  0. 8-1.0.  The  SOFC  system 
dynamic  model  is  given  in  Fig.  1. 

4.1.  Instantaneous  power  issue 

The  introduction  of  a  micro- source  within  an  existing 
power  system  assumes  that  there  may  be  a  cluster  of  small 
generators  located  in  electrical  proximity,  interconnected 
through  on-site  feeders  and  operation  can  be  isolated  or  con¬ 
nected  to  a  distribution  grid. 
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Fig.  1.  SOFC  system  dynamic  model. 


SOFCs  have  a  slow  response  to  changes  in  commands 
and  also  do  not  provide  any  kind  of  internal  form  of  energy 
storage.  These  inertia-less  systems  are  not  well  suited  to 
handle  step  changes  in  the  requested  output  power.  It  must 
be  remembered  that  the  current  power  systems  have  storage 
in  the  mechanical  energy  of  the  inertia  of  the  generators. 
When  a  new  load  is  applied  the  initial  energy  balance  is 
satisfied  by  the  inertia  of  the  system,  this  results  in  a  slight 
reduction  in  system  frequency. 

The  step  sized  power  demand  in  the  distribution  system 
should  be  instantaneously  matched  by  an  identical  supply 
of  power  from  the  micro- source  power  source. 

4.2.  Load  coming  on-line  with  grid  connection 

SOFCs  have  a  problem  in  instantaneous  power  tracking. 
SOFC  shows  a  poor  load  tracking,  as  a  load  is  applied  and 
the  micro- source  ramps  up  to  pick  up  the  whole  quota  of 
extra  power  request.  The  missing  transient  power  that  the 
micro- source  is  not  fast  enough  to  provide  is  taken  from  the 
connection  with  the  grid  that  supplies  only  a  part  of  the  local 
requested  power  during  steady- state.  Voltage  is  maintained 
by  the  power  grid. 

4.3.  Load  coming  on-line  without  grid  connection 

The  requested  power  from  the  load  coming  on-line  is  a 
step  function,  while  the  inertia-less  micro-source  always 


takes  a  finite  amount  of  time  to  ramp  up  to  the  newly 
requested  value.  Since  the  unit  cannot  change  its  output 
power  instantaneously,  the  power  is  balanced  by  voltage 
reduction.  As  the  power  injected  from  the  micro- source 
increases  to  supply  the  needed  load  the  voltage  is 
restored. 

There  is  a  missing  quota  of  energy  from  the  difference 
of  the  power  that  the  load  is  requesting  and  the  power  that 
the  unit  is  able  to  provide.  If  the  connection  to  the  main 
grid  is  active,  then  the  missing  quota  can  be  provided  by  the 
utility.  It  will  be  seen  from  the  grid  terminals  as  a  temporary, 
pulse-like  power  request. 


5.  Utility-connected  inverter  control 

For  a  six-pulse  voltage  source  inverter,  the  inverter  output 
voltage  space  vector  can  take  any  of  seven  positions  in  the 
plane  specified  by  the  ^-^-coordinates.  The  time  integral  of 
the  inverter  output  voltage  space  vector  is  called  the  “inverter 
flux  vector”  for  short.  The  flux  vector  does  not  have  the  same 
significance  as  in  motor  applications.  Rather,  it  is  a  fictitious 
quantity  related  to  the  volt-seconds  in  the  filter  inductor  [23]. 
The  d-  and  q- axes  components  of  the  inverter  flux  vector  \jrv 
are  defined  as: 

lAdv  =  f  Vd  dr  (4) 

J  —  oo 
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-f 

J  — oo 


Vqdx 


fdq  = 

The  magnitude  of  i/fv  is 

I'Avi = f 


(5) 


r2 

dv 


(6) 


The  angle  of  t/'V  with  respect  to  the  g-axis  is: 

-fdw 


5V  =  tan 


-l 


(7) 


qy 


The  d-  and  q- axes  components  of  the  ac  system  voltage 
flux  vector  t/re,  its  magnitude,  and  angle  are  defined  in  a 
similar  manner.  The  angle  between  and  V'e  is  defined  as: 

5p  =  5y  —  5e  (8) 

Control  of  the  flux  vector  has  been  shown  to  have  good 
dynamic  and  steady- state  performance.  It  also  provides  a 
convenient  means  to  define  the  power  angle  since  the  inverter 
voltage  vector  switches  position  in  the  d-q- plane,  whereas 
there  is  no  discontinuity  in  the  inverter  flux  vector. 

It  is  useful  to  develop  the  power  transfer  relationships  in 
terms  of  the  flux  vectors.  The  basic  real  power  transfer  rela¬ 
tionship  for  the  control  system  of  Fig.  2  in  the  d-q  reference 
frame  is: 


P  —  9  (< *qiq  T"  Cd^d) 


(9) 


In  (9),  eq  and  ed  are  the  q-  and  d- axes  components,  respec¬ 
tively,  of  the  ac  system  voltage  vector  E.  In  addition,  iq  and 


id  are  the  components  of  the  current  vector  /.  When  iq  and 
id  are  expressed  in  terms  of  the  fluxes,  the  equation  is  ex¬ 
pressed  as: 

3 

P  —  l^qilkqy  tyqe)  T  £j(VC/v  VC/e)]  (10) 

Taking  into  account  the  spatial  relationships  between  the 
two  flux  vectors  and  assuming  the  ac  system  voltage  to  be 
sinusoidal,  (10)  can  be  expressed  as: 

P  =  anlfe\lfy  sin  8V  (11) 

In  this  expression,  t/re  and  t/'V,  are  the  magnitudes  of  the 
ac  system  and  the  inverter  flux  vectors,  respectively,  and  5p 
the  spatial  angle  between  the  two  flux  vectors,  co  is  the  fre¬ 
quency  of  rotation  of  the  two  flux  vectors.  The  expression 
for  reactive  power  transfer  can  be  derived  in  a  similar  man¬ 
ner.  This  is 

3  co 


Q  = 


2  Lj 


COS($p  -  x//Q\ 


(12) 


Eqs.  (11)  and  (12)  indicate  that  P  can  be  controlled  by  con¬ 
trolling  5p,  which  can  be  defined  as  the  power  angle,  and 
Q  can  be  controlled  by  controlling  t/rv.  The  cross  coupling 
between  the  control  of  P  and  Q  is  also  apparent  from  these 
equations. 

The  two  variables  that  are  controlled  directly  by  the  in¬ 
verter  are  \[rv  and  <$p.  The  vector  t/rv  is  controlled  to  have  a 
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Fig.  2.  Control  system  for  the  inverter. 
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specified  magnitude  and  a  specified  position  relative  to  the 
ac  system  flux  vector  t/re.  This  control  forms  the  innermost 
control  loop  and  is  very  fast.  It  is  noted  that  both  the  in¬ 
verter  and  the  ac  system  voltage  space  vectors  are  obtained 
by  measuring  instantaneous  voltage  values  that  are  available 
locally.  The  set  points  for  the  controller  are  P*  and  E*,  and 
the  set  points  for  the  innermost  control  loop  x/r*  and  are 
derived  from  these. 

The  actual  values  of  P  and  E  calculated  from  the  feed¬ 
back  are  compared  with  the  set  values.  The  error  drives  a 
proportional  integral  (PI)  regulator,  which  generates  the  set 
points  x/r*  and  <5p  for  the  innermost  control  loop. 

The  errors  between  actual  and  desired  amounts  activate 
the  remainder  of  the  firing  scheme  only  if  they  exceed  a 
threshold  value.  If  the  error  is  larger  than  the  hysteresis  band 
(whose  widths  are  5*  and  t/r*)  then  a  decision  towards  a  new 
switching  sequence  is  made.  If  the  errors  are  within  their 
hysteresis  band,  the  switches  will  hold  their  current  status. 

The  new  switching  position  is  chosen  according  to  the 
current  angular  position  of  the  vector  at  the  inverter  termi¬ 
nals  defined  by  5V,  and  depending  on  the  fact  that  we  need 
to  increase/decrease  the  flux  magnitude  or  its  angle.  Angle 
changes  have  priority  on  flux  magnitude  changes.  Once  the 
hysteresis  threshold  is  crossed  and  the  direction  of  the  next 
switching  is  selected,  a  lookup  table  will  give  the  gating  se¬ 
quence  for  the  switches. 

Therefore,  the  SOFC  plant  has  two  major  control  loops: 

1.  Power  control:  done  by  adjusting  the  set  point  P*  of  the 
inverter  for  fast  transient  variations  and  fuel  flow  input 
control  for  slow  variations. 

2.  Voltage  control:  done  by  adjusting  the  set  point  P*  of 
the  inverter,  which  effects  the  magnitude  of  the  converter 
output  voltage. 

Using  the  fast  response  of  the  inverter,  the  maximum 
power  can  be  used  to  alleviate  transients  when  fault  occurs. 
However,  the  inverters  for  fuel  cells  have  the  following  lim¬ 
itations: 


•  Overused  fuel,  U  >  90%  (fuel  starvation  and  permanent 
damages  to  cells). 

•  Underused  fuel,  U  <  70%  (the  cell  voltage  would  rise 
rapidly). 

•  Undervoltage,  stack  voltage  <  certain  point  (loss  of  syn¬ 
chronism  with  the  network). 


Where,  the  fuel  utilization  (U)  is  the  ratio  between  the  fuel 
flow  that  reacts  and  the  input  fuel  flow,  such  as  the  follows: 


<&! 


(13) 


6.  Results 

The  analysis  is  tested  on  a  system  consists  of  13  buses 
and  is  representative  of  a  medium-sized  industrial  plant.  The 


Fig.  3.  Test  system. 

system  is  shown  in  Fig.  3  and  described  by  the  data  in 
Table  1. 

The  SOFC  plant  is  connected  at  bus  4  and  consists  of 
20  SOFCs,  as  depicted  in  Fig.  3.  The  majority  of  data  for 
the  fuel  cell  model  has  been  extracted  from  [24,25],  and  a 
commercial  leaflet  describing  a  SOFC  lOOkW  plant.  Each 
SOFC  has  the  same  parameters,  yielding  a  rated  power  of 
2  MW. 

Suppose  at  a  certain  time,  the  total  load  in  this  distribution 
system  is  as  shown  in  Table  1.  The  SOFC  plant  mainly 
provide  some  peak  shaving  capability  and  ancillary  services 
for  the  feeder. 

It  is  assumed  bilateral  contracts  to  perform  load-following 
among  SOFCs  at  a  given  feeder  is  not  allowed.  That  is,  the 
ancillary  service  must  be  coordinated  at  the  substation  or 
transmission  level. 

It  is  desirable  to  minimize  the  number  of  dynamic  equa¬ 
tions  to  be  solved  in  the  simulations.  Here,  all  SOFCs  in  the 


Table  1 


Generation,  load,  and  bus  voltage  data,  base  case 


Bus 

knag 

(%) 

P  gen 

(kW) 

Qgen 

(kvar) 

Poad 

(kW) 

Cload 

(kvar) 

1 

100.00 

7450 

540 

— 

— 

2 

99.89 

— 

— 

— 

— 

3 

99.59 

— 

— 

2240 

2000 

4 

99.64 

2000 

1910 

— 

— 

5 

101.88 

— 

— 

600 

530 

6 

99.56 

— 

— 

— 

— 

7 

101.94 

— 

— 

1150 

290 

8 

104.06 

— 

— 

1310 

1130 

9 

99.56 

— 

— 

— 

— 

10 

99.52 

— 

— 

— 

— 

11 

99.31 

— 

— 

370 

330 

12 

103.30 

— 

— 

2800 

2500 

13 

101.64 

— 

— 

810 

800 
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plant  are  regarded  as  coherent  devices,  which  have  the  same 
dynamic  response  and  share  the  generation  equally. 

6.1.  Step  change  response 


lowed  by  a  corresponding  change  in  the  fuel  flow-rate.  If  the 
fuel  cell  power  increases  while  the  fuel  flow-rate  is  main¬ 
tained  constant,  the  steady- state  cell  utilization  will  increase 
and  the  cell  voltage  will  drop. 


A  power  measurement  is  fed  in  the  fuel  flow  regulator, 
whose  output  determines  the  position  of  a  fuel  valve.  A  step 
change  in  the  input  reactant  element  will  not  be  noticed  as 
a  sudden  increase  in  charges.  Typically,  there  is  a  sudden 
but  contained  rise  in  output  power  that  takes  place  in  1-3  s, 
while  the  newly  desired  output  power  level  is  reached  lim¬ 
ited  until  the  SOFC  reaches  its  thermal  equilibrium,  which 
may  tens  of  seconds.  This  characteristic  makes  the  SOFC 
a  poor  candidate  for  isolated  systems,  where  loads  require 
instantaneous  power. 

Assume  a  SOFC  is  operating  with  constant  rated  voltage 
and  power  demand  (Pl)  0.6  p.u.  There  is  0.3  p.u.  of  step  in¬ 
crease  in  the  total  load  at  t  =  10  s.  In  the  first  2  s  or  3  s  after 
Pl  is  increased,  P  has  a  rapid  increase  due  to  the  fast  electri¬ 
cal  response  time  in  the  fuel  cell.  Subsequently,  P  increases 
slowly  and  continuously  until  reaching  the  required  power. 
This  is  due  to  the  slow  chemical  response  time  in  the  fuel 
processor.  The  total  response  time  of  P  from  0.6  to  0.9  p.u. 
is  about  30  s.  Fig.  4  shows  the  response  of  the  SOFC  to  a 
step  change  in  the  fuel  valve. 

Normally  the  fuel  cell  supports  base  load  thus,  the  fuel 
flow  regulator  in  Fig.  1  is  set  to  provide  a  constant  input 
flow-rate  for  the  cell  fuel.  The  response  of  the  fuel  flow  reg¬ 
ulator  is  slow.  The  plant,  however,  can  respond  continuously 
to  small  instantaneous  power  adjustment  by  the  control  of 
the  angle  <5*.  Any  steady-state  power  change  must  be  fol- 


6.2.  Transient  stability 


Penetration  means  the  proportion  of  the  distribution  load 
being  supplied  by  SOFC  plant  [26].  In  this  model,  an  initial 
load  of  Pi  is  assumed  and  the  penetration  is  thus, 


penetration  = 


P 

P  +  Pi 


(14) 


for  this  reason,  the  penetration  level  of  this  SOFC  plant  is 
set  at  20%. 

A  faults  occurs  that  disconnects  the  utility  supply  from 
the  distribution  system.  The  fault  clearing  is  after  0.4  s. 
The  comparison  between  Figs.  5  and  6  show  the  advan¬ 
tage  of  a  fast  power  controller.  Since  the  SOFC  inverter 
has  the  capability  to  adjust  its  angle  5*  quickly,  the  in¬ 
verter  maintains  constant  power  output  under  fast  transient 
disturbances.  In  addition,  SOFC  inverter  has  power  mod¬ 
ulation  capabilities  that  enhance  transient  stability  of  the 
system. 


6.3.  Voltage  stability 


This  section  discusses  the  nature  of  the  response  of  sys¬ 
tem  to  a  linear  increase  in  load  demand.  In  the  load  flow 
calculation,  SOFC  plant  is  considered  as  a  PV  bus  where  P 
and  E  (shown  in  Fig.  2)  are  specified. 


Fig.  4.  SOFC  plant  response  to  a  step  change  in  the  power  demand. 
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The  power  load  Pl  increases  by  a  factor  until  the  volt¬ 
age  collapse  is  reached.  The  analysis  is  tested  on  the  dis¬ 
tribution  system.  Continuation  method  is  used  to  identify 
the  location  of  bifurcation  point  (A  =  2.45),  as  presented  in 
Fig.  7. 


Fig.  8  displays  the  output  power  for  the  SOFC.  Sim¬ 
ulation  times  are  180  s  for  “mid-term”  voltage  stability 
and  1000,  10,000  s  for  “longer-term”  voltage  stability,  us¬ 
ing  a  ramp  function  equal  to  the  linear  increase  in  load 
demand. 


Fig.  6.  Fast  transient  disturbances  with  inverter. 
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0  0.5  1  1.5  2  2.5 

Load  Factor  (p.u.) 

Fig.  7.  Voltage  profiles  for  the  IEEE  13  bus  test  system. 

The  plant  can  contribute  to  the  utility  reactive  power  However,  according  to  IEEE  PI 547,  a  SOFC  system 

supply  by  continuously  adjusting  the  set  points  i/r*  °f  the  should  trip  offline  if  the  magnitude  of  the  local  electric 

inverter  to  regulate  the  magnitude  of  the  interface  bus  power  system  voltage  falls  below  88%  for  2  s  or  50%  for 

voltage.  10  cycles. 


0  20  40  60  80  100  120  140 

Time  (s) 

Fig.  8.  Dynamic  response  of  SOFC  plant.  Voltage  instability  at  180,  1000,  and  10,000  s. 
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7.  Conclusions 

The  introduction  of  SOFC  within  the  context  of  an  al¬ 
ready  existing  grid  can  provide  control  of  local  bus  voltage, 
control  of  base  power  flow  thus  reducing  the  power  demand 
from  the  main  grid  feeder,  and  ultimately  frequency  con¬ 
trol  associated  with  load  sharing  within  units  located  in  the 
micro-grid. 

The  SOFC  is  capable  of  providing  effective  load-following 
service  in  the  distribution  system.  However,  the  results  also 
show  that  the  SOFC  system  is  not  an  uninterruptible  power 
supply  and  does  not  protect  the  load  from  voltage  stability 
while  in  grid-connect  mode. 

When  SOFC  plant  is  connected  to  a  point  where  it  gives 
support  to  a  load  in  fault  conditions,  transient  stability  is 
enhanced  with  aid  of  the  SOFC  inverter. 
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